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I)  GENERAL 

This  is  a  final  report  of  contract  number  N00014-76-C-0089  on  the 
Statistical  Mechanics  of  Non-Equilibrium  Processes  in  Plasmas  and  Gases. 

II)  DESCRIPTION  OF  RESEARCH 

Since  fluids,  both  plasma  and  neutral,  are  almost  always  turbulent,  the 
turbulent  fluid  plays  an  important  role  in  a  wide  range  of  phenomena. 

The  work  in  neutral  fluids  had  been  mainly  from  the  standpoint  of  the 

fluid  equations  (Navier-Stokes  equation).  We  have  approached  the  problem  of 

obtaining  the  basic  equations  from  the  standpoint  of  statistical  mechanics* ’ ^ 

(Liouvllle  equation)  as  the  solid  foundation  of  fluid  mechanics.3)  The 

central  problem  common  to  all  approaches  is  the  closure  of  an  Infinite  set  of 

coupled  equations.  A  recent  approach^)  using  the  fluid  equations  employes  a 

three-point  Green's  function  G^3)(t,t').  The  basic  assumption  in  this  work  is 

that  the  ensemble  average  <G^3^(t,t')>  approaches  zero  as  (t-t')  >  x  where 

c 

xc  is  a  time  on  the  order  of  the  decay  time  of  the  total  energy.  This 
property  Is  used  to  obtain  closure.  A  result  of  this  theory  is  that  the 
evolution  of  the  energy  is  essentially  independent  of  the  ternary  correlation 
function  <T>  =  <v(rlt OvCr^, t)v(r.j, t)>.  We  show3)  (Appendix  A)  that  this 
conclusion  is  not  valid  in  general.  The  problem  most  likely  rests  on  the  fact 
that  <G(3)(t  t')>  is  not  independent  of  the  initial  value  of  <T>  as  was 
implicitly  assumed  in4). 

Recent  work  on  the  Adiabatic  Modifications  to  Plasma  Turbulence 
Theories^  Introduces  a  special  perturbation  theory  to  account  for  the  effects 
of  the  adiabatic  response  of  a  plasma.  We  have  shown^)  that  the  conventional 
turbulent  perturbation  theory*^,  if  properly  used,  produces  the  modified 
results  (Appendix  B). 
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III)  PERSONNEL 

The  chief  investigator  on  this  project  was  Professor  Marvin  B.  Lewis  of 
the  Department  of  Mechanical  and  Nuclear  Engineering  and  the  Department  of 
Physics  and  Astronomy  of  Northwestern  University. 

IV)  CONTRACT  PUBLICATIONS 

"Kinetic  Theory  of  Thrbulent  Flows,”  Ref.  (1) 

"Kinetic  Theory  of  Turbulent  Flows  II,”  Ref.  (2) 

"Three-Point  Method  in  Turbulence."  Ref.  (5)  (submitted  for  publication 
Physics  of  Fluids) 

Comment  on  "Adiabatic  Modifications  to  Plasma  Turbulence  Theories,” 

Ref.  (7)  (to  be  published  June  1983,  Physics  of  Fluids) 
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Three- Point  Method  in  Turbulence 


M.  B.  LEWIS 

Northwestern  University,  Evanston,  111.  (60201) 

The  three-point  aethod  for  the  description  of  turbulent  flows  rests  on 
the  assumption  that  the  average  of  the  three-point  Green's  function  approaches 
zero.  The  argument  is  made  that  the  evolution  of  the  energy  is  essentially 
Independent  of  the  Initial  value  of  the  ternary  correlation  function.  It  is 
shown  that  there  exist  Initial* conditions  for  which  this  conclusion  does  not 
hold  for  times  up  to  the  viscous  decay  time. 

PACS  47.104* 
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I.  Introduction 


Recent  work  on  the  theory  of  turbulence ^  employes  the  use  of  a  three- 
point  Green's  function  (t,t*).  The  assumption  is  made  that  the  ensemble 

average  <G^)(t,t')>  approaches  zero  for  (t-t')>Tc.  This  property,  due  to  the 
nonlinear  effects,  introduces  further  damping  in  addition  to  the  viscous 
damping  into  the  system.  On  the  basis  of  this  assumption  a  closed  equation 
for  the  energy  <V>  is  obtained.  This  equation  explicitly  contains  the  initial 
value  of  the  ternary  correlation  function  <T>  which  i6  said  to  be  arbitrary. 
The  argument  is  made  that,  since  the  dependence  on  the  initial  <T>  goes  to 
zero  due  to  the  assumed  property  of  <G^^>,  the  initial  <T>  can  be  neglected 
even  if  it  is  not  zero.  This  conclusion  is  that  the  evolution  of  <U>  is 
essentially  independent  of  the  initial  <T>. 

We  show  that  there  exist  ensembles  for  which  the  conclusion  does  not  hold 
for  a  time  tQ  from  the  initial  time  where  2t0  is  any  time  less  than  the 
viscous  decay  time. 


II.  Summary  of  the  Three-Point  Method^ 

The  Navier-Stokes  equation  for  an  incompressible  fluid  can  be  written  as 

<lt  + Vs>+l<i>t)1i<£-t)  •°  «•> 

where 

Lq  5  -uV2  ,  (lb) 

L(r,t)v(r,t)  =  /  dr  A(r,r  ):  v(r  ,t)v(r  ,t)  .  (lc) 

The  equation  for  the  energy  (covariance  correlation  function)  <U(rj,r2»t)>  = 
<v(r j,t)v(r2»t)>  for  homogeneous  systems  with  <v(r,t)>  ■  0  is 


^+Lo<E)+Lo<^o),<U<£*£o’t)>  “  -(1+tr)/dSaA(£o’Ea)!  ^E.’Xa’E’ l)>  •  <2) 
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where  tr  is  the  transpose  and  T(r , ,r0,r,,t)  b  v(r  ,t)v(r-,t)v(r,,t) ,  <  >  is 
ensemble  average  and  <T>  is  the  ternary  correlation  function.  The  basic 
problem  is  to  obtain  an  approximate  expression  for  <T>  in  terms  of  <U>  so  that 
Cq.  (2)  is  a  closed  equation. 

The  development  of  the  three-point  method  introduces  the  Green's  function 


G  defined  by 


v(r,t)  -  /dr’G(r,t;r't')*v(r\t’)  . 


Another  approach  is  to  use  the  propagator^  V  defined  by 

v(r,t)  -  W(r,t,t' )v(r,t' )  ,  (3b) 

where  W(r, t',t')  -  1.  The  equation  for  V  is  the  same  as  the  equation  for  v; 

(4r  +  L , +L)W  -  0  .  (4) 

ot  o 

The  evolution  of  T(r ^,r2,£3,t)  =  T(x,t)  is  given  by 

T(x,t)  -  W^(x,t,t'  )*T(x,t'  )  ,  (5) 

Where  W^)f  the  three-point  propagator,  is  defined  by 

W(3)(i,t,t')  i  WCtj.t.t')  W(r2,t,f)  «(t3,t,t')  ,  (6) 

where  it  is  understood  that  W(r^,t,t')  acts  only  on  functions  of  r^.  The 
equation  for  W^3^  follows  from  (4),  i.e. , 


where 


t-|^+  L3(x)  +  L3(x,t)]W<3)<x,t,t'>  -  0  , 


L3(x)  -  L  (r.)  +  L  (r,)  +  L  (r.)  , 

o  o  ~1  o  ~2  o  ~3 


L3(x,t)  -  L(r j,t)  +  L( r 2 , t )  +  L(r3,t) 


What  is  needed  in  (2)  is  <T>,  i.e.. 


<T(x,t)>  -  <VVJ'(x,t,t')T(x,t'»  .  (8) 

The  method  proceeds  with  an  expansion  of  G^33  in  terms  of  <G^33>  [Ref.  1, 


Eq.  23].  Here  we  work  with  W' 


The  corresponding  expansion  for  W'J;  is 
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V(3)(t,t')  =  <W(3)(t,t')>  +  «W(3)(t,t')  ,  (9) 

We  then  have  from  (5)  and  (9) 

T(t)  -  <W<3)(t,t’)>T(t')  +  6W(3)(t,t')T(t’)  .  (10) 

The  ensemble  average  of  T  Is 

<T(t)>  -  TIv(t)  +  Thc<t>  .  (Ha) 

TJv(t)  =  <W(3)(t,t')><T(t,)>  ,  (lib) 

T.  (t)  -  <6W(3)(t,t,)T(t,)>  .  (11c) 

nc 

The  tern  Thc  is  broken  up  into  two  terms 

Thc(t)  5  Tc(t).+  Th(t)  *  (12a) 

where  Tc  is  given  by 

Tc(t)  -  -/*,  dt’<W(3)(t,t“)><L3(t")T(t")>  ,  (12b) 


III.  The  Closed  Energy  Equation^ 

The  expression  for  <T>  has  three  parts 

<T(t)>  -  TIy(t)  +  Tc(t)  +  Th(t)  ,  (13) 

where  Tjv  is  called  the  initial  value  term  because  of  its  dependence  on  the 
initial  value  of  <T>,  Tc  is  the  cascade  term  and  Th  is  the  remainder  term. 

The  right-hand  side  of  (2)  has  corresponding  parts  SIv,  Sc,  and  Sh« 

The  assumption  employed  is  that  <W^(t2,tj)>  approaches  zero  as  (t2“tj) 
increases  with  a  characteristic  time  that  is  small  compared  with  the 
viscous  decay  time.  The  consequences  of  this  assumption  are: 

a)  The  T^  term  can  be  neglected. 

b)  That  <v  v  v  v>  that  occurs  in  S£  can  be  written  as  a  sum  of  terms  of 
the  form  <U><U>.  This  result  reduces  Eq.  (2),  except  for  the  term 
<V<3)>  in  Tc  (Eq.  12b],  to  an  equation  for  <U>.  The  closure  is 
completed  with  an  approximate  equation  for  <W^>  that  is  closed  (Ref. 
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1,  Eq.  70].  Although  no  real  estimate  was  made  for  t  it  would  appear 

c 

that  if  is  greater  than  the  decay  time  for  the  total  energy  then, 
as  in  the  usual  quasi-normal  theory"^,  the  energy  spectrum  can  become 
negative^. 

c)  The  initial  value  term  Tjy  is  negligible  after  a  time  interval  (t- 
t')>Tc.  Since  eq.  (2)  for  <U>  contains  the  initial  value  of  <T>  in 
the  term  Sjy,  it  can  be  arbitrary. 

On  the  basis  of  (c)  the  argument  is  then  made  that  TIv  can  be  neglected 
even  if  the  initial  value  of  <T>  i6  not  zero  since  it  decays  rapidly  to 
zero.  This  conclusion  is  that  the  evolution  of  <U>  is  essentially  independent 
of  the  initial  <T>.  We  shall  show  that  this  conclusion  is  not  true  in  general 
for  a  time  tQ  from  the  initial  time  where  2tQ  is  any  time  less  than  the 
viscous  decay  time.  We  compare  two  systems  that  are  related  by  velocity 
reversal.  These  two  systems  initially  have  the  same  <U>  but  one  <T>  is  the 
negative  of  the  other. 

IV.  Reversibility  of  the  Inviscld  Equations 

We  consider  the  evolution  of  a  system  from  the  initial  time  t'-  0  to  a 

time  (2tQ).  For  large  Reynolds  number  and  for  <U>  restricted  to  macroscopic 

scales  at  the  Initial  time,  the  viscous  decay  time  is  larger  than  the  decay 

time  for  the  total  energy  and  over  the  interval  (0,  2tQ)  the  effects  of 

viscosity  can  be  neglected.  Over  this  time  Interval  the  Navier-Stokes 

equation  is  reversible.  Consider  the  system  to  have  a  velocity  v+  (£.0)  at 

the  initial  time  and  a  velocity  v,  (t)  for  t  <  2t  .  Consider  a  second  system 

-«*r  O  ' 

starting  at  time  t_  with  an  initial  velocity  v  (t  )  =  -y.  (t  ).  By 
reversibility  the  evolution  of  this  second  system  over  the  interval  (t0,  2tQ) 
is  related  to  the  velocity  of  the  first  system  over  the  interval  (0,  tQ)  by 
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v  (t)  -  -v  (2t  -t)  ;  t  <  t  <  2t  (14) 

M’  O  O  O 

Calling  vn  ■  vv  -  -  we  have 

vn(t)  -  (-l)n  v°(2t  -t)  ;  t  <  t  <  2t  (15) 

-  +  o  o  o 

Now  consider  an  ensemble  of  systems,  each  member  having  a  different  v^O).  As 

in  the  three-point  method  we  consider  cases  for  which  <$v+(t)>  -  0.  At  tc 

construct  a  second  ensemble  formed  by  taking  for  each  member  of  the  ensemble 

v  (t  )  ■  — v  (t  ).  It  follows  from  (15)  that 
—  o  ~+  o 


<vn(t)>  -  (-1)°  <y°(2t  -t)>  (16) 

t  <  t  <  2t 
o  o 

Special  cases  of  (16)  for  n*2  and  3  are 

<D  (t)>  -  <0  (2t  -t)>  ;  <T  (t)>  -  -  <T  (2t  -t)>  (17) 

—  +  o  —  +  o 

We  compare  the  evolution  of  the  two  ensembles  starting  from  tD  over  the 
interval  (tQ,  2tQ).  At  the  initial  time  t  **  te,  from  (17),  both  ensembles 
have  the  same  energy  <U_(t0)>  “  <U+(t0)>  and  the  ternary  correlation  functions 
are  the  negative  of  each  other  <T_(tQ)>  ■  -<T+(t0)>.  From  (16)  with  t  -  tQ 
all  even  moments  are  equal  and  all  odd  moments  are  the  negative  of  each 
other.  Assuming  that  the  ensemble  at  tQ  does  not  weight  v  and  -v  equally 
(e.g.  Gaussian)  these  two  ensembles,  having  different  initial  moments,  evolve 
in  different  ways.  The  evolution  of  the  second  ensemble  over  (tc,  2tQ)  is 
related  to  the  evolution  of  the  first  ensemble  over  (0,  tc)  by  (16).  If  the 
first  ensemble  at  tQ  were  e.g.  Gaussian  then  velocity  reversible  doe6  not 
produce  a  new  ensemble  i.e.,  the  first  and  second  are  the  same  ensemble.  From 
the  moment  point  of  view  all  odd  moments  are  zero  and  therefore  all  moments  of 
the  two  ensembles  are  equal.  Assuming  that  this  is  not  the  case  the  evolution 
of  the  two  ensembles,  even  though  they  start  with  the  same  <U>,  can  evolve  in 
distinctly  different  ways.  Suppose  that  the  first  ensemble  evolves  over 


•  M 
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(0,  2t0)  In  such  a  way  chat  the  energy  flows  from  large  to  small  scales  and 
therefore  from  large  to  small  scales  over  (tQ,  2tQ).  The  second  ensemble 
energy  over  the  Internal  (te,  2tQ),  by  (17),  flows  from  small  to  large 
scales.  The  usual  inviscid  quasi-normal  theory^  explicitly  shows  this  type 
of  behavior.  The  above  argument  is  general  in  that  it  does  not  Involve  a 
specific  theory. 

For  the  three-point  theory,  <w£3^(t,tQ)>  and  <W^3^(t,tQ)>  go  to  zero  for 

t  ■  t  +  t  .  According  to  (a)  and  (c)  the  right-hand  side  of  the  energy 
o  c 

equation  (2)  for  t  >  tQ  +  xc  depends  only  on  T£  and  according  to  (b)  Tc  is  a 
functional  of  <U>.  For  the  <U>  of  the  two  ensembles  to  evolve  from  tQ  in  the 
same  way  T^Ct)  must  essentially  equal  Tc_(t).  For  the  first  ensemble  T^  is 
given  by 

-  -/*  dt"<W<3)(t,t“)XL3(t-)T+(t“)>  ;  to  <  t  <  2to  (18) 

o 

The  contribution  to  the  Integral  comes  from  an  interval  (t—rc)  <  t"  <  t  i.e. , 
it  depends  on  <U+(t“)>  for  t"  in  the  neighborhood  of  t"  -  t.  For  the  second 
ensemble  Tc_(t)  is 

T  (t)  -  -L  dT<W(3)(t,T)XL3(T)T  (t)>  ;  t  <  t  <  2t  (19) 

c—  t  —  -  —  o  o 

o 

3 

By  reversibility  and  the  fact  that  L  is  linear  in  v  we  have 
<L3(t)T  ( t)>  -  <L3 (2t  -t )T . (2t  -t )> 

—  —  TO  TO 

so  that 

T  (t)  -  -  /°  dt"<W(3)(t,2t  -t")XL3(t")T .  (t“)>  (20) 

C~  2t  -t  0  +  + 

o 

The  contribution  to  the  integral  comes  from  an  interval  (2t^-t)  <  t"  < 
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(2tQ-t)  +  tc  i.e.,  it  depends  on  <U+(t")>  for  t"  in  the  neighborhood  t"*2t0~t. 
For  a  system  (the  first  ensemble)  that  is  evolving  over  an  interval  (0,  2tQ) 
the  value  of  <U+>  for  a  time  in  the  interval  (0,  tQ)  is  different  from  its 
value  for  a  time  in  the  Interval  (t0,  2tQ)  and  therefore  T^  is  in  general 
different  than  Tc_. 

Our  conclusion  is  that  for  a  given  Initial  <U>  there  do  exist  ensembles 
(related  by  velocity  reversible)  of  differing  initial  <T>  that  evolve  in 
distinctly  different  ways  over  a  time  interval  tQ.  The  above  argument 
restricts  the  initial  <T>  for  which  the  three-point  method  is  valid.  There 
can  be  a  class  of  ensembles  with  the  same  initial  <U>  but  different  initial 
<T>  that  evolve  in  essentially  the  same  way  over  the  Interval  tQ.  However  the 
class  cannot  contain  pairs  of  ensembles  that  are  related  by  velocity  reversal. 
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Comment  on 


"Adiabatic  Modifications  to  Plasma  Turbulence  Theories" 

M.  B.  Lewis 

Northwestern  University,  Evanston,  Illinois  60201 

In  Ref.  1  Catto  introduced  a  new  perturbation  scheme  to  account  for  the 

adiabatic  response  of  a  plasma  and  obtained  the  nonlinear  dielectric  function. 

Krommes^  calculated  the  dielectric  function  from  the  conventional  theory 

2  2 

obtaining  a  result  differing  from  Catto's  by  a  term  (k  Ur  )  and  then  remarked 
that  "Catto  makes  subsldary  approximations  which  lead  him  to  neglect  the 
(k^Dr  )  term”.  Catto3'  responded  that  the  term  (k*Dr  )  appears  erroneously  in 
Ref.  2  "because  the  properties  of  the  adiabatic  response  are  not  properly 
preserved  by  the  average  trajectories".  The  purpose  of  this  note  is  to  show 
that  Catto's  new  scheme  can  be  obtained  from  the  conventional  scheme  if  use  is 
made  of  the  constraint 

lfm  -  0  (1) 

where  FM  is  the  part  of  the  distribution  function  that  depends  on  the  energy, 

L"L+L' ,  L£ v*  | — ,  L'5  —  E'»  and  E'“V#  is  the  fluctuation  in  the  field.  The 
~  3r  m  ~  3v  ~  T 

ensemble  average  field  E  is  zero. 

We  start  with  Weinstock' s'  formulation  of  the  conventional  scheme  for  f, 

the  fluctuating  part  of  the  distribution  function  F 

t  _ 

f(r,v,t)  -  U.(t,0)f(0)  -  /  dt'UA(t,t* )L' (t' )F(t' )  (2) 

AS  M  A  Q  A 

where  F-F+f.  The  operator  is  defined  by 

llt  +  FlVt-t.>‘0  «> 

where  P=-(A-1)L  and  A  is  an  operator  which  takes  the  ensemble  average  of 
everything  on  which  it  operates.  Equation  (2)  is  exact  and  therefore  can  be 
employed  in  the  problem  of  the  near  equilibrium  plasma.  In  this  case^ 


-1- 


-2- 


F-Fm+G ,  FM»Foexp['|i] ,  F0-c  exp[-mv2/2T] .  We  also  have  f-F^+g  where 
Fj^fF^-F^)  and  g ~(G-G).  Equation  (2)  la  then 

f  -  UA(t,0)[F^(0)+g(0)]  -  /  dt'UA(t,t')L’G  -  I 


I  =  /  dt’U^t.t'jL'FJt') 

0  A  M 

The  approximation  developed  by  Weinstock^'  for  UA  Is 

0A(t,t’)(l-A)  -  (l-A)U(t.t')  (5) 

where  U  is  the  ensemble  average  of  the  Vlasov  operator  U.  The  problem 
referred  to  by  Catto^  Is  in  using  U  in  1  since  FM  satisfies  (1).  This 
problem  can  be  circumvented  by  a  transformation  on  I.  The  integral  I  can  be 
written  as 

f  —1  — 

I  -  J  dt'UA(t,0)UA  (t*  ,0)L'  (t*  )Fjj(t*  ) 

—1  ® 

where  UA  is  the  inverse  of  UA  and  satisfies  the  equation 

5TuI1<t'to)  -  "*  ('•to)p  *  0  (6> 

From  rFM-0  it  follows  that  L  F^—L'F^  and  from  this  that  L'FM«- PF^.  Using 
this  relation  and  (6)  in  I  we  have 

t  3F'(t') 

I  -  UA(t,0)F^(0)  -  F^(t)  +  /  dt'UA(t,t’>  — —  (7) 

Equation  (4)  is  then,  using  (5) 

t  _  3F*(t’) 

f  -  U(t,0)gC0)  +  F'(t)  -  /  dt’U(t,t' )[— TP — +  1>'C]  (8) 

“  0  __ 
Equation  (8)  is  equivalent  to  Catto's^  eq.  (6b)  except  that  the  operator  U 

3F' 

acts  on  the  entire  term  -—T  and  not  just  the  |*-r  part.  Using  the  same 
approximation  as  Catto,  following  his  eq.  (6b)  we  have*7 


?k  s  /  d2fk  -  f  ♦k(t)U+i“/dvdvoUk  jU(v,vo)Fo(vo)]  - 


*  ¥  ♦  d20F0^20^1+  /  dt(iu))expli((i)-kvaT)-k2Dr  3/3j} 


where  U.  is  defined  in  ref.  (2).  The  conventional  result  is  obtained  from 
k,w 

(2)  using  (S)  and  replacing  F  by  FQ  and  results2^  in  (9)  with  the  factor 


2  2 

(la)  replaced  by  (iut-k  Dt  ).  The  conventional  perturbation  (2)  in  conjunction 
with  the  conatraint  (1)  gives  the  result  of  Catto  (8)  and  produces  in  the 
ezplicit  fora  for  f  a  term  differing  by  (k  Dt  )  from  the  explicit  form  for  f 
obtained  only  from  the  conventional  perturbation  (2). 
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